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ABSTRACT 

Galaxy clusters exhibit regular scaling relations among their bulk properties. These relations estab- 
lish vital links between halo mass and cluster observables. Precision cosmology studies that depend 
on these links benefit from a better understanding of scatter in the mass-observable scaling relations. 
Here we study the role of merger processes in introducing scatter into the M-Tx relation, using a 
sample of 121 galaxy clusters simulated with radiative cooling and supernova feedback, along with 
three statistics previously proposed to measure X-ray surface brightness substructure. These are the 
centroid variation (w), the axial ratio (77), and the power ratios (P20 and P30). We find that in this 
set of simulated clusters, each substructure measure is correlated with a cluster's departures <51nTx 
and S In M from the mean M-Tx relation, both for emission- weighted temperatures Tew and for 
spectroscopic-like temperatures TgL, in the sense that clusters with more substructure tend to be 
cooler at a given halo mass. In all cases, a three-parameter fit to the M-Tx relation that includes 
substructure information has less scatter than a two-parameter fit to the basic M-Tx relation. 
Subject headings: galaxies: clusters: general, X-rays: galaxies: clusters 



1. INTRODUCTION 

Clusters of galaxies play a critical role in our under- 
standing of the Universe and its history and are po- 
tentially powerful tools for conducting precision cosmol- 
ogy. For example, large cluster surveys can discrimi- 
nate between cosmological models with different dark- 
energy equations of state by providing complementary 
observations of the shape of the cluster mass function, 
evolution in the number density of clusters with red- 
shift, and bias in the spatial distribution of clusters 
(IWang fc Steinhardtl 1199a iLeyine et all 120021: |Hul l2003t 
iMaiumdar fc Mohrll2004j : fVoitll2005f l. However, this po- 
tential to put tight constraints on the properties of dark 
energy will be realized only if we can accurately measure 
the masses of clusters and can precisely characterize the 
scatter in our mass measurements. 

Scatter in X-ray cluster properties is directly related to 
substructure in the intracluster medium. If clusters were 
all structurally similar, then there would be a one-to-one 
relationship between halo mass and any given observable 
property. Generally speaking, deviations from a mean 
mass-observable relationship are attributed to structural 
differences among clusters. One kind of structural differ- 
ence is the presence or absence of a cool core, in which 
the central cooling time is less than the Hubble time at 
the cluster's redshift, and the prominence of a cool core is 
observed to be a source of scatter in sc a ling relationships 
tiFabian et all 1 1994 iMarkevitchl [19971 IVoit et all I2002 : 
iMcCarthv et al.l l2004f) . We also expect structural dif- 
ferences to arise from substructure in the dark matter, 
galaxy, and gas distributions. For instance, there may be 
a spread in halo concentration at a given mass, variations 
in the incidence of gas clumps, differences in the level of 
AGN feedback, or various effects due to mergers. All of 
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these deviations can be considered forms of substructure 
that produce scatter in the mass-observable relations one 
would like to use for cosmological purposes. While it may 
ultimately be possible to constrain the amount of scat- 
ter and its evolution with re dshift using self-calibration 
techniques ([Lima fc Hull2005t ), such constraints would be 
improved by prior knowledge about the relationship be- 
tween scatter and substructure. 

Traditionally, the most worrisome form of substructure 
has been that due to the effects of merger events. Clus- 
ters are often identified as "relaxed" or "unrelaxed" , with 
the former assumed to be nearly in hydrostatic equilib- 
rium and the latter suspected of being far from equilib- 
rium. Cosmological simulations of clusters indicate that 
the truth is somewhere in between. The cluster popula- 
tion as a whole appears to follow well-defined virial rela- 
tions with log-normal scatter around the mean, showing 
that clusters do not cl eanly separat e into relaxed and 
unrelaxed systems (e.g. lEvrard et ail (|2007t ) ). Even the 
most relaxed-loo king clusters are not qu ite in hydrostatic 
equilibrium (e.g. lKravtsov et all (|2006l )). Instead of sim- 
ply being "relaxed" or "unrelaxed," clusters occupy a 
continuum of relaxation levels determined by their re- 
cent mass-accretion history. 

Quantifying this continuum of relaxation offers oppor- 
tunities for reducing scatter in the mass-observable rela- 
tions. If mergers are indeed responsible for much of the 
observed scatter around a given scaling relation, then 
there may be correlations between a cluster's morphol- 
ogy and its degree of deviation from the mean relation. 
Once one identifies a morphological parameter that cor- 
relates with the degree of deviation, one can construct a 
new mass-observable relation with less scatter by includ- 
ing the morphological parameter in the relation. Such 
an approach would be analogous to the improvement of 
Type la supernovae as distance indicators by using light- 
curve stretch as a sec ond parameter to indicate the su- 
pernova 's luminosity (|PhillipsHl99"llRiess et alll 1996ft . 

Here we investigate how merger-related substructure 
in simulated clusters affects the relationship between a 
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simulated cluster's mass and the t emperature of its in - 
trac luster medium build i ng upon IBuote fc Tsail (|1995l ) 
and lO'Hara et alj (|2006l ). IBuote & Tsail (|1995l) auanti- 
fied the morphologies and dynamical states of observed 
clusters and found structure measures to be an indica- 
tor of the dynamical state of a cluster. lO'Hara et alJ 
(2006) also examined morphological measurements, for 
both observed and simulated clusters, and found that 
simulations without cooling showed no correlation be- 
tween substructure and scaling relation scatter. In this 
work we examine substructure for simulated clusters with 
radiative cooling and focus on the idea that merger pro- 
cesses introduce intrinsic scatter into the M-Tx relation- 
ship by displacing clusters in the M-Tx plane away from 
the mean X-ray temperature (Tx)| M at a given mass 
M, either to higher or lower avera ge ICM temperature . 
We then adopt a s et of statistics (jBuote fc Tsail Il995t 
lO'Hara eT al. 2006) for quantifying galaxy cluster sub- 
structure and merger activity in order to investigate this 
hypothesis. Section 2 discusses the M-Tx scaling rela- 
tionship in our sample of simulated clusters and shows 
that disrupted-looking clusters in this sample tend to 
be cooler at a given cluster mass. In Section 3 we at- 
tempt to quantify the relationship between morphology 
and temperature using four different substructure statis- 
tics and compare it to similar studies. We then show 
that substructure in these simulated clusters indeed cor- 
relates with scatter in the M-Tx relationship and assess 
the prospects for using that correlation to reduce scatter 
in the M-Tx plane. Section 4 summarizes our results. 

2. MASS-TEMPERATURE RELATION IN SIMULATED 
CLUSTERS 

This study is based on an analysis of 121 clus- 
ters simulated using the cosmo logical hydrod ynamics 
TREE+SPH code GADGET-2 (|Springe]|l2005t ). which 
were simulated in a standard A cold dark matter 
(ACDM) universe with matter density Qm = 0.3, h = 
0.7, fif, = 0.04, and erg = 0.8. The simulation treats ra- 
diative cooling with an optically-thin gas of primordial 
composition, includes a time-dependent UV background 
from a population of quasars, and handles star forma- 
tion and supernova feedback using a two-phase fluid 
model with cold star-forming clouds embedded in a hot 
medium. All but f our of the clusters ar e from the simu- 
lation described in Bor gani et al.l (|2004h . who simulated 
a box 192 h^ 1 Mpc on a side, with 480 3 dark matter par- 
ticles and an equal number of gas particles. The present 
analysis considers the 117 most massive clusters within 
this box at z = 0, which all have M200 greater than 
5 x 10 13 h~ 1 M ( r ) . By convention, Ma refers to the mass 
contained in a sphere which has a mean density of A 
times the critical density p c , and whose radius is denoted 
byi? A 

That cluster set covers the ~ 1.5-5 keV temperature 
range, but the 192 /i -1 Mpc box is too small to contain 
significantly hotter clusters. We therefore supplemented 
it with four clusters with masses > 10 15 hr 1 Mq and tem- 
peratures > 5 keV drawn from a dar k-matter-only sim - 
ulation in a larger 479 hr 1 Mpc box (|Saro et al.ll2006f ). 
The cosmology for this simulation also was ACDM, but 
with erg = 0.9. These were then re-simulated including 
hydrodynamics, radiative cooling, and star formation, 
again with GADGET-2 and using the zoomed-initial- 



conditions technique of iTormenl (|1997t ). with a fourfold 
increase in resolution. This is comparable to the reso- 
lution of the clusters in the smaller box. Adding these 
four massive clusters to our sample gives a total of 121 
clusters with M200 in the interval 5 x 10 13 ft. _1 M© to 
2 x 10 15 h~ l M Q . 

We first need to specify our definitions for mass and 
temperature. In this paper, cluster mass refers to A/2oo- 
For temperature, we use two definitions. The first is the 
emission-weighted temperature 

/T[n 2 A(T)]ri 3 x 
TeW " J^A(T)d^x ' (1) 

where n is the electron number density and A(T) 
is the usual cooling function for intracluster plasma. 
The second is the s pectroscopic-like temperature of 
iMazzotta etail (|2004h . 

_ fT[n 2 T-V 4 } d 3 x 

T SL- J n 2 T -3/4 d 3 x > ( 2 ) 

where the power-law weighting function replacing A(T) 
is chosen so that TgL approximates as closely as possible 
the temperature that would be determined from fitting 
a single-temperature plasma emission model to the inte- 
grated spectrum of the intracluster medium. The pres- 
ence of metals in the ICM of real clusters introduces line 
emission that co mplicates the co mputation of TgL for 
clusters <3 keV (jVikhlininl 2006). However, the simu- 
lated spectra for the clusters in our sample are modelled 
with zero metallicity, which eases this restriction in our 
analysis. 

Figure [T] shows the mass-temperature relations based 
on these definitions for our sample of simulated clusters. 
The best fits to the power-law form 

have the coefficients Mq ~ 5.9 x 10 13 h~ 1 M Q ,a ~ 
1.5 for Tx corresponding to Tsl and Mq ~ 4.4 x 
10 13 h~ 1 M Q ,a ~ 1.4 for T x corresponding to Tew As is 
generally the case for simulated clusters, the power-law 
indices of the mass-temperature relations found here are 
consis tent with cluster self-similarity a nd the virial the- 
orem ([Kaiserlll986t iNavarro et al.lll995f ). These relation- 
ships have scatter, which we characterize by the standard 
deviation in log space <7i n M about the best-fit mass at 
fixed temperature Tx- When relating M to the emission- 
weighted temperature Tew 5 we find cr\ n M — 0.102. 
When relating cluster mass M to the spectroscopic-like 
temperature TgL, the scatter is o\ n M = 0.127. That the 
scatter is larger for the spectroscopic-like temperature is 
not surprising, given the sensitivity of Tsl to the thermal 
complexity of the ICM. 

Figure Q] also highlights two subsamples for each defi- 
nition of temperature, selected based on the clusters' de- 
viations in In Tx space from the mean mass-temperature 
relation. In each panel, open circles represent the eight 
clusters that have the largest positive deviations and are 
therefore "hotter" than other clusters of the same mass, 
while filled circles represent the eight with the largest 
negative deviation and are "cooler" than other clusters of 
the same mass. In general, these temperature estimates 
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Fig. 1. — Mass-temperature (M-Tx) relationships for the 121 
clusters in our sample. Spectroscopic-like temperature Tsl is plot- 
ted in the top panel, while emission-weighted temperature Tew 
is plotted in the bottom panel. The solid lines show the best-fit 
power-law relation M = Mo(Tx/3 keV) a over the whole sample, 
while the dashed lines show the best-fit relation for systems with 
Tx > 2 keV. Open circles represent the clusters that have the 
greatest positive temperature offset from the mean relationship, 
and filled circles represent the clusters with the greatest negative 
temperature offset. 

are well correlated, so that hotter outliers in Tew 
also hotter outliers in Tsl and cooler outliers in Tew are 
also cooler outliers in Tsl- Since Figure Q] distinguishes 
the most extreme outliers for the two temperature esti- 
mates, this distinction may define slightly different sets, 
though they still overlap. 

Figure [2] presents a gallery of surface brightness maps 
for two sets of eight clusters with the most extreme off- 
sets from the mean M-Tsl relation. The eight unusually 
hot clusters are in the top panel, and the eight cooler 
clusters are in the bottom panel. In these plots the hot- 
ter clusters appear more symmetric, and are seemingly 
"more relaxed," and the cooler clusters appear less sym- 
metric and seemingly "less relaxed." The gallery as a 
whole therefore suggests that relaxed clusters tend to be 
hot for their mass and unrelaxed clusters tend to be cool 
for their mass. 

At first glance, the result that disrupted-looking clus- 
ters in cosmological simulations tend to be cooler than 
other clusters of the same mass may seem counterintu- 
itive, since one might expect that mergers ought to pro- 
duce shocks that raise the mean temperature of the in- 
tracluster medium. This fin ding has also been n o ted b y 
iMathiesen fc Evrardl (|2001h and iKravtsov et alj (|2006l ). 
Cluster systems in the process of merging tend to be cool 
for their total halo mass because much of the kinetic en- 
ergy of the merger has not yet been therma l ized. 

The idealized simulations of iPoole et alJ (|2007f ) illus- 
trate what may happen to the ICM temperature during 
a single merger. Before the cores of the two merging sys- 
tems collide, the mean temperature is cool for the over- 



FlG. 2. — Surface-brightness contour maps for sixteen of the clus- 
ters in our sample, overlaid with a circular aperture of radius -R500 ■ 
The top panel has eight maps representing the clusters that have 
the largest positive deviation 8 InTx from the mean M-Tx relation, 
calculated using the spectroscopic-like temperature TgL- These are 
the objects represented by open circles in the top panel of Figure 
[T] The bottom panel has eight maps representing the clusters that 
have the largest negative deviation from the same M-Tx relation- 
ship, represented by the filled circles in Figure [T] Note that the 
"hotter" clusters in the top panel appear more symmetric, while 
the "cooler" clusters in the bottom panel appear more irregular. 

all halo mass because it is still approximately equal to 
the pre-merger temperature of the two individual merg- 
ing halos. There is a brief upward spike in temperature 
when the cores of the merging halos collide, after which 
the system is again cool for its mass. Then, as the re- 
maining kinetic energy of the merger thermalizes over 
a period of a few billion years, the temperature gradu- 
ally rises to its equilibrium value. The merging system 
therefore spends a considerably longer time at relatively 
cool values of mean temperature for its halo mass than 
at relatively hot values. Hence, such simulations sug- 
gest a possible explanation for why more relaxed sys- 
tems would tend to lie on the hot side of the M-Tx 
relation, while disrupted systems would tend to lie on 
the cool side. A caveat, however, is that the current 
generation of hydrodynamic cluster simulations tend to 
produce relaxed clusters whose temperature profiles con- 
tinue to rise to smaller radii t han is observed in real 
clusters (jTornatore et all l2003t iNagai et al.ll2007t ). po- 
tentially enhancing average temperatures for such sys- 
tems. As a separate test of this effect, we excise the core 
regions from our sample clusters, calculate new substruc- 
ture measures and new emission-weighted temperatures 
for the core-excised clusters, and repeat our analysis. 

3. QUANTIFYING SUBSTRUCTURE 

The question we would like to address in this study 
is whether the surface-brightness substructure evident 
in Figure [2] is well enough correlated with deviations 
from the mean mass-temperature relation to yield use- 
ful corrections to that relation. In order to answer that 
question, we need to quantify the surface-brightness sub- 
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structure in each cluster image, so that we can deter- 
mine the degree of c orrelation across the entire sample. 
lO'Hara et al.l (|2006f ) explored the relationship between 
cluster structure and X-ray scaling relations in both ob- 
served and simulated clusters, and we adopt their suite 
of substructure measures in this study. These include 
centroid varia t ion, a xial ratio, and the power ratios of 
iBuote fc Tsail (|1995f ). In this section we define and dis- 
cuss those statistics and apply them to surface-brightness 
maps made from three orthogonal projections of each 
cluster. Then we assess how well these statistics correlate 
with offsets from the mean mass-temperature relation. 

3.1. Axial Ratio 

The axial ratio r\ for a cluster surface-brightness map is 
a measure of its elongation, which is of interest because 
it has been found from simulations that t he ICM is of- 
ten highly elongated dur ing merger events (|Evrard et al.l 
l!993t|Pearce et all 1 9941 ). It is computed from the second 
moments of the surface brightness, 

Mij = ^IxXiXj. (4) 

The summation is conducted over the coordinates 
(x\, X2) of the pixels that lie within an aperture centered 
at the origin of the coordin ate system to which {x\ 1 X2) 
refer. Following the work of lO'Hara et alJ ()2006l ). we use 
an aperture of radius -R500 centered on the brightness 
peak. We then compute 77 from the ratio of the non-zero 
elements that result from diagonalizing the matrix M. 
That is, 

D = U T MU, (5) 
where U is a diagonalizing matrix for M, and 

HI: « 

The axial ratio is therefore defined to be in the range 
77 G [0,1], with 77 = 1 for a circular cluster. Of course 
there are other choices for the origin of the coordinate 
system, besides using the brightness peak. For instance, 
in order to avoid misplaced apertures yielding artifi- 
cially low axial ratios for nearly circular distributions, 
one could adjust the position of the aperture to seek a 
maximum in 77. Doing this, we sometimes find that 77 « 1 
even for non-circular clusters, as is evident in Figure [3l 
This figure depicts the surface-brightness map of what 
appears to be a disturbed cluster, chosen from among 
those in our sample that appear by eye to be the most 
unrelaxed. Yet, it happens to have an axial ratio very 
close to 1 for an aperture placed so as to maximize 77. 
This example demonstrates that, while the axial ratio 
statistic may yield results consistent with a visual in- 
terpretation of cluster substructure, it is also capable of 
unexpected results for some clusters. 

To further illustrate this point, we have computed an 
axial ratio value for this cluster for every possible choice 
of aperture placement. Apertures of radius i?soo were 
centered on each and every pixel within the surface- 
brightness map, provided the aperture so placed does 
not reach the edge of the map. This procedure gener- 
ated an axial-ratio "surface" mapping all of the aperture 
placements to a value of 77. Figure [4] shows the axial-ratio 




Fig. 3. — Surface-brightness contour plot of an asymmetric clus- 
ter which, for certain choices of aperture placement, yields an axial 
ratio close to 1. The circle represents an aperture of R500 




Fig. 4. — Surface of axial ratio r] as a two-dimensional function 
of the coordinates of the aperture center. The axial ratio statistic 
appears to be ill-defined for this cluster. 
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Fig. 5. — Abstract surface of axial ratio 77 as a two-dimensional 
function of the coordinates of the aperture center. This is a relaxed 
cluster, for which the axial ratio is better-defined. 

surface for the cluster in Figure [31 For comparison pur- 
poses, Figure [5] presents an axial-ratio surface map for 
a very symmetric, uniform, and apparently relaxed clus- 
ter, in which the cluster's brightness peak reassuringly 
corresponds to the aperture location that maximizes 77. 
In contrast, the presence of two peaks in the axial-ratio 
surface for the asymmetric cluster shows that 77 can some- 
times depend strongly on aperture placement. Ideally, we 
would like to place the aperture on the "center" of this 
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Fig. 6. — Relationship between a cluster's deviation <51nTx from the mean M-Tx relationship and four measures of substructure: the 
centroid variation w, the axial ratio r], and the power ratios P20 and P30. The lower panels are for spectroscopic-like temperature, and 
the upper ones are for emission-weighted temperature. The light gray bands and vertical dashed lines show the extremes in substructure 
measurements, with 80% of our clusters falling in the central region between the bands. The black filled circles correspond to clusters above 
2 keV, while the plus signs correspond to clusters below 2 keV. Finally, the solid lines indicate the best-fitting linear relationships between 
our substructure measures and temperature offset, for the above 2 keV sub-sample denoted by the black filled circles. 



cluster, but the center of an unrelaxed cluster can be 
difficult to define, meaning that the axial ratio statistic 
may be likewise ill-defined for such clusters. 

3.2. Power Ratio 
The power-rati o statistics ()Buote fc Tsail Il995t 



lO'Hara et al.ll2006h quantify substructure by decompos- 
ing the surface-brightness image into a two-dimensional 
multipole expansion, the terms of which are calculated 
from the moments of the image, computed within an 
aperture of radius R ap : 



i (.Rap) 



£(x')(R') r 



cos mo 



fx' 



(7) 



b m (R. 



R'<B. a 



£(£' )(i?') m sin m0' d 2 x'. (8) 



The power in terms of order m is then 

0*m + &m) 
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For m = 0, the power is given by 

P = [aoMRap)] 2 . (10) 
The power ratios P m o = P m /Po are then dimensionlcss 
measures of substructure which have differing interpreta- 
tions. For instance, P10 quantifies the degree of balance 
about some origin and can be used to find the image cen- 
troid, P20 is related to the ellipticity of the image, and 
P30 is related to the triangularity of the photon distri- 
bution. As in the case of the axial ratio computations, 
we set the aperture radius P ap equal to Psoo- The most 
appropriate place to center the aperture is at the set of 
pixel coordinates that minimizes P10, which we achieve 
using a self-annealing algorithm. 



3.3. Centroid Variation 

The centroid variation statistic w is a measure of the 
center shift, or "skewness", of a two-dimensional pho- 
ton distribution. It is measured for a cluster surface- 
brightness map in the following way. For a set of surface- 
brightness levels one finds the centroids of the corre- 
sponding isophotal contours and computes the variance 
in the coordinates of those centroids, scaled to i?soo- Here 
we select 10 isophotes evenly spaced in log/x between 
the minimum and maximum of Ix within an aperture 
of radius P500 centered on the brightness peak, so as 
to adapt to the full dynamic range of surface bright- 
ness for different clusters. We employed this adaptive 
scheme because using one set of isophotes for all clusters 
tended to ignore important substructure in less massive 
clusters when they had surface brightness substructure 
inside P500 but outside of the lowest isophotc. 

3.4. Substructure and Scaling Relationships 

Using the quantitative measures of substructure de- 
scribed in the previous section, we can test the signifi- 
cance of the relationship between substructure and tem- 
perature offset hinted at in Figure [2] We begin by treat- 
ing four of our substructure statistics — centroid variation 
w, axial ratio 77, and power ratios P20 and P30 — as dif- 
ferent imperfect measurements of an intrinsic degree of 
substructure S. Figure [5] shows the relationship between 
substructure and a cluster's deviation SlriTx from the 
mean M-Tx relation for each substructure measure. In 
each case we present results for both the spectroscopic- 
like temperature Tsl and the emission-weighted temper- 
ature Tew- Note that centroid variation w and the power 
ratios P20 and P30 have large dynamic ranges, whereas 
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the axial ratio rj is always of order unity. We therefore 
attempt to fit the relationships between JlnTx and the 
different substructure measures with the following forms: 



SlnTx = 



Aw a 
B + (3r) 
CP20 

np x 



(11) 



To visually indicate where the bulk of our substruc- 
ture measures lie, Figure LH has light gray bands cov- 
ering the extremes, so that 80% of our sample clus- 
ters have substructure measures lying between the ex- 
tremes. The power ratios in our study generally span 
two decades (in units of 10 -7 ), from ~2 — 300 for F20 
and from ~0.01 — 1 for P30. Th e se ran g es are consis- 
tent w ith t hose of iBuote fc Tsai| (|1995l ). lO'Hara et all 
<|2006f >. and lJeltema et alJ (|2007f ). The measurements of 
axial ratio in our sample, with 80% of clusters having 
77^0.4 — 0.95, cover a slightly wide r rang e than do the 
simulated clusters of lO'Hara et al.l (|2006l ). Finally, our 
measurements of centroid variation, with 80% of clusters 
having w[i?5q n ]~0.0 1 — 0.1, are again similar to those of 
lO'Hara et all (j200l . 

As denoted in Figure [6] by black filled circles, the sys- 
tems with Tx above 2 keV occupy a slightly narrower 
range of substructure values than the systems below 2 
keV, which are denoted by plus signs. For the axial ra- 
tio and the power ratios, the variance is 15 to 25 per- 
cent larger among the low-temperature systems when 
compared to the systems with Tx > 2 keV. For cen- 
troid variation the variance among the low-temperature 
systems is approximately the same as it is among the 
high-temperature systems. However, it is not clear that 
there is a significant correlation between substructure 
and mass, since the mean substructure values are gener- 
ally very similar between the low-temperature and high- 
temperature subsamples. The mean value of the power 
ratio P30 is significantly larger for the low-temperature 
subsample, however this measure also has the weakest 
correlation with offsets from the mean M-Tx relation. 

To test whether the low-mass clusters in our sample 
significantly boost the overall scatter in the M-Tx rela- 
tion, we perform a cut at 2 keV and fit this relation both 
to the whole sample and to the sub-sample above 2 keV. 
Figure [7] shows the residuals in mass, actual minus pre- 
dicted, where the predicted mass derives only from the 
M-Tx relation. The plus signs indicate clusters whose 
mass is predicted from an M-Tx relation derived from 
all 121 clusters. The black filled circles indicate clusters 
that are above 2 keV in X-ray temperature, with the mass 
estimated using the sub-sample M-Tx relation. There is 
a negligible reduction in scatter, from 0.127 to 0.124 for 
Tsl and from 0.102 to 0.094 for Tew, suggesting that at 
best only a modest improvement is found in our sample if 
we remove the low-mass systems. In order to test the de- 
gree to which incorporating substructure measures adds 
to this modest improvement, when we compare mass es- 
timates derived using substructure to those derived only 
from the M-Tx relation, we focus on clusters above 2 
keV in the rest of our analysis. 

Figure LH] shows that for our simulated clusters, a 
greater amount of measured substructure tends to be as- 
sociated with "cooler" clusters while less substructure 
tends to be associated with "hotter" clusters. Also, the 
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Fig. 7. — Comparison of mass offset 5\nM between true mass 
and predicted mass, based on the Af(Tx) relation. Plus signs indi- 
cate residuals for masses estimated from the M-Tx relation derived 
from all 121 clusters, while black filled circles are for masses esti- 
mated from the M-Tx relation for clusters above 2 keV. Upper 
panels are for spectroscopic-like temperature and lower panels are 
for emission-weighted temperature. The standard deviations a in 
the residuals are given in each plot. 

centroid variations w are more highly-correlated with 
SlnTx than are the other substructure parameters. We 
interpret this to mean that the centroid variation is a 
better predictor of the offset in the M-Tx relationship 
than are the power ratios and the axial ratio, though all 
four measures appear to be related to the temperature 
offset. Again, in this figure we denote systems above 2 
keV by black filled circles, and systems below 2 keV by 
plus signs. 

Correlations between substructure and S In Tx can po- 
tentially be exploited to improve on mass estimates of 
real clusters derived from the mass-temperature rela- 
tion. Instead of computing the temperature offset at 
fixed mass, we can determine a substructure-dependent 
mass offset at fixed temperature and then apply it as a 
correction to the predicted mass M py:et x(Tx) one would 
derive from the mean M-Tx relation alone. To assess 
the prospects for such a correction, based on this sample 
of simulated clusters, we first define the mass offset from 
the mean mass-temperature relation to be 



5\nM(T x ) = In 



M 



M prcd (Tx) 



(12) 



where M is the cluster's actual mass, and examine the 
correlations between substructure measures and ShiM. 
Figure [8] shows the results. These plots show mass pre- 
dictions from both the M-Tsl relation and the M-Xew 
relation. Consistent with our analysis of (SlnTx, the cen- 
troid variation w appears to be a more effective predictor 
of the mass offset S In M(Tx). Nonetheless, all four mea- 
sures of substructure appear to be correlated with mass 
offset. 
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Fig. 8. — Relationship between substructure and mass offset S In M(Tx) from the mean M-Tx relationship for the same substructure 
measures as in Figure [6] The lower panels are for Tx = TgL, and the upper panels are for Tx = Tew As in Figure [6] the solid lines indicate 
the best-fitting linear relationships to the above 2 keV sub-sample denoted by black filled circles, and the gray bands and vertical dashed 
lines mark the extremes in substructure between which 80% of our clusters lie. Also as in Figure [6] the plus signs correspond to systems 
below 2 keV., while the black filled circles correspond to systems above 2 keV. 



In order to incorporate a substructure correction into 
the mass-temperature relation, we perform a multiple 
regression, fitting our simulated clusters' mass, temper- 
ature, and substructure data to the form, 

logM pred (T x , S) = logMo + odogTx + (3S, (13) 
where S represents one of the following substructure 
measures: 77, logiu, logf^o, or log-Pso- This fit gives us a 
substructure-corrected mass prediction M pre d(Tx, S) for 
each substructure measure, and we can assess the effec- 
tiveness of that correction by measuring the dispersion 
of the substructure-corrected mass offset 

M 



5 In M(T X , S) = In 



LM pred (T x ,S) 



(14) 



between the revised prediction and the true cluster mass. 

Figure [9] shows the results of that test. Open circles in 
each panel indicate mass offsets <Hn M(Tx) without sub- 
structure corrections, which have a standard deviation 
cm(t)- Filled circles indicate mass offsets S In M(Tx, S) 
with substructure corrections, which have a standard de- 
viation a M (t.s) ■ The upper set of panels shows results for 
Tsl, and the lower set is for Tew- I n each case, incorpo- 
rating a substructure correction to the mass-temperature 
relation reduces the scatter, yielding more accurate mass 
estimates. The centroid variation corrections are the 
most effective, reducing the scatter in mass from 0.124 
to 0.085 in the M-T SL relation and from 0.094 to 0.072 
in the M-Tew relation, though admittedly this is again 
a modest improvement. Although non-negligible struc- 
ture correlates significantly with offsets in the M-Tx 
plane, apparently it does so with substantial scatter. 
This scatter may be partly due to projection effects, in 
which line-of-sight mergers are discounted by the mea- 
sures of substructure and may dilute their corrective 
power peltema et al.ll2007l ). 



Lastly, Figure [TU] shows the results for a similar analy- 
sis to that of Figure [9l except that in this case we have 
excised a region of radius 0.15i?5oo around the center of 
each cluster and recomputed Tew- We do this to test 
whether offset in temperature, whose correlation with 
substructure is the basis of our correction scheme, stems 
from a potentially unrealistic feature, which is that the 
cores of many real clusters have temperature profiles that 
decline at larger radii than occurs in simulated clusters. 
As in Figure [H we restrict our analysis to clusters above 
2 keV. After doing this test, for Tew excising the core ac- 
tually increases the scatter in M-Tx from 0.094 to 0.106. 
It may be that by removing the bright central region, the 
average temperature becomes more sensitive to structure 
outside the core. Also, this figure shows that the effect of 
incorporating substructure measurements into the mass- 
estimates is still present. The scatter is reduced to 0.075 
for to, 0.093 for r), 0.090 for P 20 and 0.094 for P 30 . Figure 
[TOl summarizes the results of this test, which support the 
conclusion that the reduction in scatter we realize using 
substructure is a real effect and not an artifact of known 
defects in the simulations. 

3.5. Comparisons with Other Substructure Studies 

O Ti ara et al.l (|2006l ) examined the relationship be- 
tween galaxy cluster substructure and X-ray scaling re- 
lationships, including the M-Tx relation, using both a 
flux-limited sample of nearby clusters and a sample of 
simulated clusters, and found a greater amount scatter 
among the more relaxed clusters in their observed sam- 
ple. Contrasting that result they also found a greater 
amount of scatter among the more disrupted clusters in 
their simulation sample, though they characterize the ev- 
idence for this second result to be weak. Finally, they see 
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Fig. 9. — Comparison of mass offset <51nAf between true mass and predicted mass, based on the M(Tx) relation (open circles) and the 
M(Tx, S) relation (filled circles). Upper panels are for Tsl and lower panels are for T^w- The standard deviations a in the residuals are 
given in each plot. 



no evidence in either sample for more disrupted clusters 
to be below the mean, and the more relaxed clusters to 
be above. One difference between our study and theirs 
is the presence of radiative cooling and supernova feed- 
back in the simulation that produced our cluster sample. 
Also, the focus of our work is different from theirs in that 
we concentrate on the degree of correlation between the 
amount of substructure and the size and direction of the 
offset from the mean relation. We do find significant ev- 
idence of this correlation, such that relaxed clusters are 
hotter than expected given their mass. We also test, as 
best we can given our simulation sample, the hypothesis 
that substructure can be used to improve mass estimates 
derived from the ICM X-ray temperature. It is possible 
that our detection of a correlation between substructure 
and temperature offset arises from the additional physics 
in our simulated clusters, since when radiative cooling is 
included, cool lumps may be better preserved than in 

simulations that don't include cooling. 

Our results are in agreement with IValdarnin i (2006), 
who examined substructure in clusters simulated with 
cooling and feedback and found that unrelaxed clusters, 
identified with a larger power ratio P30, have spectral- fit 
temperatures biased low relative to the mass-weighted 
temperatures. This trend aligns with our finding that 
the spectroscopic-like temperature Tsl is lower than the 
best-fit temperature at fixed mass for clusters with larger 
power ratios P 2 o aim -£30- However, iValdarninil (2006) 
did not investigate the effectiveness of substructure mea- 
sures in reducing scatter in the mass-temperature rela- 
tion. 

Our res ults are also i n agre ement with some of the 
results of Ueltema et al.l (|2007j ). who have recently in- 
vestigated correlations between substructure and offsets 
in mass predictions in simulated clusters. They found 
that measuring cluster structure is an effective way to 



correct masses estimated using the assumption of hy- 
drostatic equilibrium, which tend to be underestimates. 
Our findings support these results, given that we find 
substructure can be used to correct masses estimated di- 
rectly from the M-Tx relationship. There also are dif- 
ferences between our findings and theirs. They report 
that the M-Tx relation for their simulation sample shows 
no dependence on structure, whereas the clusters in our 
sample exhibit offsets that correlate with the degree of 
substructure. One possibility is that these differences 
stem from differences in the simulations' feedback mech- 
anisms. Another possibility is that some of the offset we 
observe derives from enhanced temperatures in simula- 
tions with radiative cooling. As we describe in section 
3, we perform a test in which we estimate Tew using 
projected surface-brightness and temperature maps, in 
order to remove the core regions from our analysis, but 
this may be less effect ive than properly exc ising the cores 
i n the simulation s . as | Jeltema et al. (2007) have done. 

iKravtsov et alj (|2006l ) also looked at the relationship 
that cluster structure has to the M-Tx relation in simu- 
lated clusters, to show that the sensitivity of mass proxies 
Yx and Ygz to substructure is not very strong. They 
divided their sample into unrelaxed and relaxed sub- 
samples, based on the presence or absence of multiple 
peaks in the surface-brightness maps of clusters, and 
found the normalization of the M-Tx relation to be bi- 
ased to cooler temperatures for the unrelaxed systems. 
Other workers also have looked at the relationship be- 
tween the M-Tx relation and subs tructure, as reflected 
in the X-ray spectral properties. iMathicscn fc Evrardl 
(|2001h have examined the ratio of X-ray spectral-fit tem- 
peratures in hard and full bandpasses for an ensemble of 
simulated clusters, and found it to be a way of quanti- 
fying the dynamical state of a cluster. We consider our 
approach of using surface-brightness morphology infor- 
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Fig. 10. — Comparison of mass offset 51nM between true mass and predicted mass, based on the M(Tx) relation (open circles) and 
the M(Tx, S) relation (filled circles). These results are for are for T^w, with the central 0.15 -R500 region removed both from the average 
temperature and from the substructure measures. The standard deviations a in the residuals are given in each plot. 

M-Tew relation was 0.094. Including substructure as a 
third parameter reduced the scatter to 0.072 for centroid 
variation, 0.084 for axial ratio, 0.081 for P20, arid 0.084 
for P30. Scatter about the basic two-parameter AT- Tsl 
relation was 0.124, and including substructure as a third 
parameter reduced the scatter to 0.085 for centroid vari- 
ation, 0.112 for axial ratio, 0.110 for P20, and 0.108 for 
P30. As one last test, and to increase our confidence 
that our substructure measures are not relying on po- 
tentially non-physical core structure in the simulations, 
we also repeated the comparison of mass-estimates for 
Tew, with the core regions of the clusters excised. First, 
removing the core slightly increased the scatter in M-Tx 
possibly by making the average temperature more sen- 
sitive to structure outside the core. Second, even with 
the cores removed the improvement in mass-estimates 
obtained using substructure information remains. Based 
on these results, it appears that centroid variation is the 
best substructure statistic to use when including a sub- 
structure correction in the M-Tew relation. However, 
the correlations we have found in this sample of simu- 
lated clusters might not hold in samples of real clusters, 
because relaxed clusters in the real universe tend to have 
cooler cores than our simulated clusters do. 



mation to be co mplementary to theirs. More recently, 
iKav et aTl (|2007f ) performed an interesting analysis on 
another large-volume simulation sample, using as sub- 
structure metrics the centroid variation and measures 
of concentration to report evolution in the luminosity- 
temperature relationship. Specifically, they report that 
the more irregular clusters in their sample lie above the 
mean M-Tx relation (i.e., they are cooler than average), 
for the spectroscopic-like temperature Tsl- 

4. SUMMARY 

Using a sample of galaxy clusters simulated with cool- 
ing and feedback, we investigated three substructure 
statistics and their correlations with temperature and 
mass offsets from mean scaling relations in the M-Tx 
plane. First, we showed that the substructure statistics 
w , Vi P20 and P30 all correlate significantly with (51nTx, 
though with non-negligible scatter. In all cases this scat- 
ter is larger for <51n Tsl than it is for <51n Tew- Next, we 
considered the possibility that M-Tx scatter is driven by 
low-mass clusters. We tested the degree to which scatter 
can be reduced by filtering out these systems. This con- 
sisted of performing a cut at 2 keV, for which we saw that 
it yielded a modest improvement in mass estimates. To 
see whether incorporating substructure could refine these 
mass estimates, we first showed that w, rj, P20, and P30 
correlate significantly with the difference S In M between 
masses predicted from the mean M{Tx) relation and 
the true cluster masses, with non-negligible scatter that 
again is less for M(Tew) than it is for M(Tsl)- Then we 
adopted a full three-parameter model, M-Tx-S, which 
includes substructure information S estimated using w, 
rj, P20, and P30. Scatter about the basic two-parameter 
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